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We report on a measurement of the top-quark electric charge in tt events in which one W boson originating 
from the top-quark pair decays into leptons and the other into hadrons. The event sample was collected by the 
CDF II detector in ^/s = 1.96 TeV proton-antiproton collisions and corresponds to 5.6 fb _1 . We find the data 
to be consistent with the standard model and exclude the existence of an exotic quark with —4/3 electric charge 
and mass of the conventional top quark at the 99% confidence level. 



PACS numbers: 14.65.Jk, 12.15Ji, 14.65.Ha 

I. INTRODUCTION 

Since the discovery of the top quark (f) fjHH], the CDF and 
DO collaborations, joined recently by the LHC experiments, 
have measured several of its properties to be consistent with 
standard model (SM) predictions. Determining that the top 
quark decays into a W + boson and a bottom quark (b), while 
the anti-top quark decays to a W~ boson and an anti-bottom 
quark would ensure indirectly that the electric charge of the 
(anti-)top quark is indeed (— )2/3 as expected in the SM. If 
events were found to contain decays into aW" and bottom- 
quark final state, the charge of the decaying particle would be 
-4/3, incompatible with the SM top quark. Motivation for a 
measurement was proposed in Ref. [3], where such a hypoth- 
esis was put forward. In this model, an exotic quark of mass 
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around 170 GeV/c 2 is assumed to be part of a fourth gen- 
eration of quarks and leptons, while the standard-model top 
quark is heavier than 230 GeV/c 2 . Even though this model 
is by now strongly disfavored by other measurements [4, 5], 
the charge correlations between jets initiated by b or b quarks 
and W bosons in tt events have not yet been definitively es- 
tablished. The existence of an exotic decay combination (b 
coupled to W~ and b coupled to W + ) has already been con- 
strained experimentally ||6l |7D, but with less sensitivity than 
the present measurement. 

In this article we analyze tt candidate events and treat the 
SM and exotic-quark hypotheses exclusively. We analyze tt 
candidate events in the final state containing hadrons from the 
decay of one W boson and an electron or muon and corre- 
sponding antineutrino from the decay of the other W boson. 
We first determine the charge of the W boson (using the charge 
of the lepton or the opposite charge for the hadronically de- 
cayed W boson). Then we pair the W boson with the jet orig- 
inating from a b quark (b jet) from the same top-quark decay. 
Finally we determine the charge of the b jet using an opti- 
mized jet-charge algorithm, JetQ ll8l-fTTll. Pairings where the 
charge of the W boson is opposite to the JetQ value are classi- 
fied as standard-model-like (SM-like) decays, while pairings 
where the charge of the W boson is of the same sign are clas- 
sified as exotic-model-like (XM-like) decays. 

In Sec. [II] we briefly describe the CDF II detector. The 
data sample and event selection are presented in Sec. [Ill] and 
Monte Carlo simulations in Sec. [IV] Section [VI discusses the 
method to pair the W boson with the correct b jet. The JetQ 
algorithm used to assign a charge to the b jet, as well as its cal- 
ibration using data, are described in Sec.[VT] The backgrounds 
and the possible biases they may induce in the measurement 
are investigated in Sec. lVIII In Sec. |VIII| the systematic uncer- 
tainties are presented while Sec. [IX]explains how the pairing 
purity and JetQ purity are combined to obtain the signal pu- 
rity, i.e., the probability of correctly identifying a signal event 
as coming from the SM or the XM. The statistical treatment 
of the data is described in Sec.[X]and the results are presented 
and discussed in Sec.lXIl 
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II. THE CDF II DETECTOR 

The CDF II detector is described in detail in Refs. ifTlfll . 
The subdetectors most relevant to this measurement are briefly 
described in this section. The detector is approximately her- 
metic over the full angular coverage and is composed of 
a charged particle tracker embedded in an axial magnetic 
field of 1.4 T, surrounded by electromagnetic and hadronic 
calorimeters and muon detectors. A cylindrical coordinate 
system with z-axis directed along the proton beam is used. 
The polar angle is defined with respect to the proton beam 
direction and (]) is the azimuthal angle about the z-axis. Pseu- 
dorapidity is defined as r| — — lntan(8/2). 

The charged particle tracker is composed of silicon micro- 
strip detectors lfl4l - [l6ll covering the pseudorapidity range of 
|r| I < 2 and providing 1 1 pm spatial resolution on each mea- 
surement point in the r-(|) plane, crucial for the identification 
of secondary vertices characteristic of jets originating from b 
quarks. The silicon detectors are surrounded by a 3.1m long 
open-cell drift chamber lfl7ll . which measures the momenta 
of charged particles within a pseudorapidity range of |r| | < 1 . 
The calorimeter covers the pseudorapidity range of |r|| < 3.6 
and is segmented into projective towers that point towards the 
nominal center of the interaction region. The electromagnetic 
portion is a lead-scintillator sampling calorimeter [18], which 
also contains proportional chambers and resistive strips at a 
depth corresponding to the typical maximum shower intensity 
for electrons. The hadronic portion is an iron-scintillator sam- 
pling calorimeter [19]. Muon detectors are located outside the 
calorimeters. Two sets of drift chambers separated by steel 
absorber, the CMU JH and CMP JUIH, cover the pseudo- 
rapidity range |r|| < 0.6, and layers of drift tubes sandwiched 
between scintillation counters, the CMX 112 U 12211 . cover the 
range 0.6 < |r|| < 1.0. 



III. DATA SAMPLE AND EVENT SELECTION 

This analysis is based on a data sample corresponding to 
an integrated luminosity of 5.6 fb _1 collected with the CDF 
II detector between February 2002 and February 2010. The 
events first have to pass an inclusive-lepton online event se- 
lection (trigger) that requires an electron with Ej > 1 8 Ge V 
or muon with pj >18 GeV/c [23]. We then select events of- 
fline with a reconstructed isolated electron Ej (or muon pr) 
greater than 20 GeV (GeV/c), and missing E t {#t) >20 GeV 
02411 . In addition we require events to have at least four jets, 
three of them with Ej >20 GeV and |r|| < 2.0 and another jet 
with Ej > 12 GeV and |r|| < 2.4. We explicitly reject events 
that have two or more leptons to ensure that the final sample 
does not include events where both W bosons decay into lep- 
tons (dilepton channel). 

The electron selection relies on the accurate geometrical 
match between a reconstructed track and some energy de- 
position in the electromagnetic calorimeter. We also require 
that the amount of energy deposited in the hadronic calorime- 
ter be significantly less than in the electromagnetic calorime- 
ter. An isolation criterion requires the transverse energy in 



the towers not assigned to the electron, within a cone of 
AR = \J (Ar|) 2 + (A(j)) 2 =0.4 centered around the lepton, to 
be less than 10% of the candidate electron Ej. 

In the muon selection, a track candidate from the tracker 
is matched to a track segment (stub) in one or more of the 
muon drift chambers. We require either a stub in both the 
CMU and CMP chambers, or a stub in the CMX chamber, 
and refer to the resulting muon candidates as CMUP or CMX 
muons respectively. The energy deposited in the region of 
the calorimeter to which the trajectory of the candidate muon 
extrapolates is required to be consistent with the expectation 
for a minimum-ionizing particle. The isolation criterion for 
muons, similar to that for electrons, is that the calorimeter 
transverse energy in a cone of AR = 0.4 around the extrapo- 
lated muon track (not including the muon energy deposition 
itself) must be less than 10% of the muon pj. Details on the 
electron and muon identification are discussed in Ref. ll25ll . 

The muon acceptance is increased by approximately 20% 
by including events containing muons that cannot be triggered 
on directly. Such events must pass a different trigger, which 
requires a missing transverse energy larger than 35 GeV and at 
least two jets of Ej > 10 GeV. Candidates are selected if they 
contain a CMX stub in a region not covered by the inclusive 
lepton trigger, or a stub only in the CMU or CMP chambers, or 
an isolated track not fiducial to any muon detector. Muons in 
these categories, called extended muons, are also required to 
pass the isolation criterion and to have pj > 20GeV/c. Dilep- 
ton veto and jet requirements are the same as those applied 
to events selected from the inclusive lepton trigger. To ensure 
full efficiency of the trigger, the extended muon candidates are 
also required to have two jets with Ej > 25 GeV, one of which 
should be central (|r|| < 0.9) and separated from the other by 
Afl>1.0. 

The jet reconstruction is based on a calorimeter-tower- 
clustering cone algorithm with a cone size of AR = 0.4. Tow- 
ers corresponding to selected electrons are removed before 
clustering. The observed Ej for jets is corrected for the ef- 
fects of jet fragmentation, calorimeter non-uniformities and 
the calorimeter absolute energy scale ll26ll . 

Due to the presence of a neutrino leaving the detector unde- 
tected, there will be an imbalance in the transverse momentum 
of the event. Consequently, events are expected to have some 
missing transverse energy $t, and we require fyj > 20 GeV. 

The data set selected above, called "lepton+jets" (LJ), is 
dominated by QCD production of W bosons with multiple jets 
("W + jets"). To improve the signal-to-background ratio we 
identify events with two or more b jets, i.e., we require at least 
two of the jets to contain a secondary vertex, characteristic of 
a B hadron having decayed. This secondary vertex algorithm 
is tuned such that the efficiency of identifying a b jet is about 
50%, and results in a probability of about 2% of misidentify- 
ing a light-quark jet. More information about this algorithm 
can be found in Ref. 1251. 
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IV. MONTE CARLO SIMULATION 

The tt Monte Carlo (MC) simulation used in this measure- 
ment relies on PYTHIA version 6.216 lf27ll for event genera- 
tion and parton showering. The top-quark mass used is 172.5 
GeV/c 2 . Samples generated with other values of the top- 
quark mass are studied for any dependence of the measure- 
ment on this parameter. A sample of tt events generated with 
HERWIG version 6.510 l28ll is used to estimate a possible sys- 
tematic uncertainty due to the choice of generators. Most of 
the background samples rely on PYTHIA except for the W + 
jets background, which is generated using transition matrix el- 
ements calculated by ALPGEN version 2.10' JH] and PYTHIA 
for parton showering. Parton distribution functions are mod- 
eled with CTEQ5L BOIl . The interactions of particles with the 
detectors are modeled using GEANT3 ||3l]|, and the G FLASH 



parametrization [32l] for showers in the calorimeters. Details 
on the implementation and tuning of the CDF II detector sim- 
ulation are found in Ref. i33"ll . 



V. PAIRING BETWEEN THE W BOSON AND THE b JET 

Each event contains a lepton, multiple b-jet candidates, and 
non-b jets. In order to assign the four highest-pr jets to the 
four final-state quarks from the tt decay and to associate the 
lepton with the b jet from the decay of the top quark that 
produced the leptonically-decaying W boson, we use the top- 
quark mass kinematic fitter described in Ref. [34], which min- 
imizes a % 2 variable that incorporates constraints on the top- 
quark mass m t , fixed at 172.5 GeV/c 2 , and on the W -boson 
mass mw, fixed at 80.42 GeV/c 2 . The % 2 is given by 



x 2 = 



E 

i=t,4jets 



i \2 



(Pt-Pt) 



(Pj ~Pj ) (mjj-mw) 1 (m (y -m w ) 2 {m hj j-m t ) 2 (m Mv -m t ) 2 



(1) 



]=x,y 



The first term evaluates the difference between the best-fit 
value (pr) and the observed value (pr) of the transverse mo- 
mentum for the four highest-/;^ jets and the lepton. The sec- 
ond term evaluates the difference between the best-fit and the 
observed value of the unclustered energy, which represents the 
energy in the calorimeter towers not associated with the jets 
or primary lepton. The last four terms represent the mass dif- 
ferences between the W boson and its decay products and be- 
tween the top quark and its decay products. The parameter m t 
is not floating, in contrast to Ref. 113411 . The a,- and a, variables 
are the uncertainties on the observed momenta values, while 
Ow represents the decay width of the W boson (2.12 GeV/c 2 ), 
and Gt is the quadrature sum of the theoretical width of the 
top quark (1.5 GeV/c 2 ) and the experimental uncertainty on 
its mass (0.9 GeV/c 2 ). Since events may contain two, three, 
or four jets identified as b jets by the secondary-vertex algo- 
rithm, there are two, six, or twelve possible assignments of b 
jets to W bosons, respectively. For each W — b pairing two % 2 
values are computed to allow for the unknown z component of 
the neutrino momentum. Choosing the combination that min- 
imizes this x 2 leads to a purity p va ir (the probability of correct 
W-b pairing) of 76%, as estimated with the PYTHIA tt MC 
sample. By imposing an upper threshold to the value of the 
minimum % 2 , the purity is increased but the event selection 
efficiency is reduced. We identify the optimal configuration 
by maximizing eD 2 obtained from the tt simulated sample, 
where e is the efficiency of the % 2 requirement and D is the 
dilution, defined as D = 2p pa i r — 1. By restricting the analy- 
sis to events in which the minimum % 2 does not exceed 9, we 
achieve an efficiency on signal of 53.2 ± 0.1% with a purity 

Ppair Of 83.3 ±0.1%. 



VI. CHARGE OF A b JET 

We use the jet-charge (JetQ) algorithm to determine which 
of the high-pr b jets characteristic of a tt event originated 
from a b quark, and which from a b quark. We select tracks 
with impact parameter M35I1 less than 0.15 cm with respect to 
the primary vertex and pj larger than 1 .5 GeV/c within a cone 
of AR < 0.4 around the b jet axis. We only compute JetQ if 
there are at least two such tracks within this cone. We then 
sum up the charges of those tracks with weights that depend 
on their momentum component along the jet axis: 



JetQ 



HiPtrack ' Pjet) ' Qtrack 



YJ\.P track ' Pjet) 



0.5 



(2) 



where pj et (p trac k) is the momentum vector of the jet (track) 
and Quack is the charge of the particle associated to the track. 
Track requirements and the choice of the 0.5 exponent result 
from an optimization of JetQ on the simulated tt sample. If the 
JetQ value is positive we assign the bottomjet to a b quark, if it 
is negative we assign the bottomjet to a b quark. Monte Carlo 
studies indicate that this algorithm has a selection efficiency 
of 97.9±0.1% and apurity per identified b jet of about 60. 8± 
0.1%. 



A. Calibration of the JetQ purity in data 

Since the simulation does not model the jet fragmentation 
reliably, we correct the purity of the JetQ algorithm obtained 
from the simulation by using a dijet data sample enriched in 



6 



20< E T (AJ) <35 GeV 
50< E T (AJ) <75 GeV 
E T (AJ) > 90 GeV 

Bottom 




0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



. 20< E T (AJ) <35 GeV 
. 50< E T (AJ) <75 GeV 
E T (AJ) > 90 GeV 

Charm 



P Trel (GeV/c) 



0.5 1 1.5 2 2.5 3 3.5 4 

P Trel (GeV/c) 



> 

0) 

O 

IT) 
CM 
O 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



20< E T (AJ) <35 GeV 

50< E,.(AJ) <75 GeV 

■ ■■■ E T (AJ)>90GeV 

Light 




0.5 1 1.5 2 2.5 3 3.5 4 



P Tre| (GeV/c) 



20< E T (AJ) <35 GeV 
50< E T (AJ) <75 GeV 
E T (AJ) > 90 GeV 

Bottom 




12 3 4 

M vtx (GeV/c 2 ) 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



. 20< E T (AJ) <35 GeV 
. 50< E T (AJ) <75 GeV 
E T (AJ) > 90 GeV 

Charm 



2 3 4 

(GeV/c 2 ) 



> 
o 
CD 



> 
LU 



0.4 




0.35 




0.3 




0.25 


: i 


0.2 




0.15 




0.1 




0.05 





20< E T (AJ) <35 GeV 

50< E,.(AJ) <75 GeV 

■ ■■■ E T (AJ)>90GeV 

Light 




2 3 4 

(GeV/c 2 ) 



FIG. 1 : Examples of pjrel (top) and secondary- vertex mass (bottom) templates for bottom, charm, and light quarks for different values of 
away -jet Ej. 



heavy flavor. This sample is collected with a trigger that re- 
quires a central muon with pj > 8 GeV/c. Events are required 
to have a tag muon jet with Ej > 20 GeV that contains a 
muon with pj > 9 GeV/c inside the cone, and a probe away 
jet with Ej > 20 GeV and with A(j) > 2 with respect to the 
muon jet. We require both jets to be identified as b jets using 
the secondary-vertex algorithm, but a more selective variant 
of the tagger is used for the muon jet. The JetQ purity is ob- 
tained as the fraction of selected events in which the charge of 
the muon is opposite to the JetQ value of the away jet. The ob- 
served purity is corrected for a number of effects. If the muon 
originates from a b —> c —> p. cascade decay, its charge is the 
opposite of the one it would have if coming directly from a b 
decay (secondary fraction); if the B meson undergoes mixing, 
the charge of the muon may also flip sign (mixing fraction); 
and finally, if one of the two b jet candidates is misidentified, 
no correlation between the JetQ value and the charge of the 
identified muon (non-M? fraction) is present. The first two 
effects can be obtained from simulation. The last effect is cal- 
culated from a fit of simulation to data. 

In order to obtain the bb fraction of the dijet sample, we 
use two independent fits. We first extract the b fraction in 
muon jets by fitting the distribution of prreh the component of 
the muon momentum transverse to the jet direction, which is 
enhanced at larger values for muons originating from b quark 
jets. FigureQ](upper panels) shows a selection of the pj re i tem- 
plates used. For this fit we combine the charm and light-quark 
templates since they are very similar. Then, we determine the 



b fraction in away jets by fitting the secondary-vertex-mass 
distribution, which is enhanced at higher values when the par- 
ent quark is heavier. Figure [T](lower panels) shows a selection 
of secondary-vertex mass templates used; the template shapes 
depend on the away -jet Ej. To allow for the possibility that 
the simulated sample might not model reliably the Ej distri- 
bution of light quarks that are misidentified as b quarks, we 
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FIG. 2: Ratio between the purity of the JetQ algorithm in data and 
Monte Carlo simulation as a function of away -jet Ej. The dotted line 
is the result of fitting a constant through the data points. 

perform all template fits in nine independent ranges of away- 
jet Ej. Since the b fractions of the muon and away jets are 
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obtained from independent fits, we have no information on 
their correlation in the dijet sample. However we can obtain 
the highest (lowest) value of the bb fraction by assuming that 
this correlation is maximal (minimal). We then estimate the 
bb fraction in each Ej range as the average of the upper and 
lower limits in the range, and set the corresponding uncer- 
tainty to half the difference between the limits. 

Combining the bb fraction with the secondary and mixing 
fractions we correct the bias in the measured purity in each 
away-jet Ej bin. We compute a scale factor SFjq as the ra- 
tio of the purity obtained in the dijet data sample to that ob- 
tained in a corresponding simulated sample. We see no de- 
pendence of the scale factor on the away-jet Ej, as shown 
in Fig. 12 We estimate a total systematic uncertainty on the 
JetQ scale factor of 3.2%, coming from uncertainties on the 
template shape (2.3%), the fit strategy (1 .8%), and the Ej de- 
pendence (1.4%). We obtain a value of the scale factor of 
SFjq = 0.99 ± 0.01 (stat) ± 0.03 (syst). 



VII. BACKGROUNDS 

In the following, signal refers to events with either a SM 
tt pair, or a pair of exotic quarks with mass 172.5GeV/c 2 . 
The exotic quarks are simulated using the standard tt Monte 
Carlo described in Sec.[lV] The dominant background is QCD 
production of W plus multijet events. These events enter the 
signal sample when two of the jets are b jets (W+HF), or light 
quark jets are misidentified as b jets (mistag). Other back- 
grounds include QCD multijet events where a jet is misidenti- 
fied as a lepton and two jets are b jets or misidentified as such 
(QCD fakes), single-top-quark events, and diboson events. 
The amount of background is moderate (« 15%) because at 
least two jets are required to be identified as b jets. 

We obtain the background predictions with the same 
method as for the cross-section measurement of Ref. lf36ll . We 
compute the efficiency of the % 2 requirement and JetQ selec- 
tion using Monte Carlo simulation for each background with 
the exception of the QCD fakes, for which we use data. Fi- 
nally, we search for correlations between the charge of the 
primary lepton and the JetQ value of the corresponding b jet 
in each background source. This correlation is expressed as 
the fraction of the total number of W-b pairs that are clas- 
sified as SM-like. We expect this fraction to be 50%, i.e., 
the same probability for pairs to be SM- or XM-like, except 
for two processes, single-top-quark production and QCD bb 
production where a lepton from the semileptonic b decay is 
misclassified as primary lepton. For the first process we rely 
on the simulation to estimate the possible correlation, while 
for the second process we use a data sample where all the LJ 
selection requirements are applied except those of the lepton 
selection, and we require instead that the lepton fail at least 
two identification criteria. The composition of this sample is 
dominated by QCD background events. Table Q] summarizes 
the signal and background predictions. Table [TTI summarizes 
the amount of correlation for each background. Background 
sources for which no effect is expected are assigned a corre- 
lation of 0.5. The signal correlation (purity) is defined in Sec. 



US 

In Fig. [3] we show the % 2 distribution used to assign the 
lepton to the correct jet, while in Figs. 2]and[5]we show the 
distributions of the number of tracks in the JetQ calculation 
and the lepton pj, respectively. 




i 2 



FIG. 3: Distribution of minimum % 2 for events passing selection re- 
quirements described in Sec. [in] Shaded histograms show signal and 
background predictions stacked to form the total prediction. The ar- 
row shows the % 2 upper threshold. 



■ W+HF 




Number of tracks in Jet Q 

FIG. 4: Distribution of number of tracks entering the JetQ calcu- 
lation. Shaded histograms show signal and background predictions 
stacked to form the total prediction. The purity of the JetQ algorithm 
is calibrated as described in Sec. |VI Al and a scale factor to the Monte 
Carlo is obtained to account for modeling discrepancies. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties come from modeling of the ge- 
ometrical and kinematic acceptance, knowledge of the sec- 
ondary vertex tagging efficiency, the effect on the acceptance 
of the uncertainty on the jet energy scale, uncertainties on the 
background predictions, and the uncertainty on the luminos- 
ity. 
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TABLE I: Background and signal expectations before and after the % 2 and JetQ criteria (columns 2 and 5). The efficiencies of these criteria 
are shown in columns 3 and 4. The last column includes a factor of two because each selected event contains two W — b pairs, each providing 
a "quark candidate" (SM- or XM-like candidate). The uncertainties are discussed in Sec. I Villi 



Process 


Events 
before criteria 


% z requirement 
efficiency 


JetQ 
efficiency 


Quark candidates 
after criteria 


W+HF 


66 ±22 


0.150±0.004 


0.970 ±0.003 


19.5±6.4 


QCD fakes 


18 ± 14 


0.17 ±0.08 


0.88±0.12 


5.4±4.8 


Diboson 


4.7±0.7 


0.22 ±0.02 


0.97 ±0.01 


2.0±0.4 


Mistag 


9.7±2.6 


0.15 ±0.02 


0.96 ±0.02 


2.8±0.8 


Single top 


10.6± 1.3 


0.210±0.004 


0.970 ±0.003 


4.4±0.5 


Total background 


109 ±26 






34±8 


Signal 


670 ±110 


„ „ 9 ±0.001(stat) 
U - 3JZ ±0.005(syst) 


Q7Q ±0.000(stat) 
U - y/V ±0.002(syst) 


700 ±120 



TABLE II: Correlation between lepton charge and JetQ in background and signal events. The last two columns show the expected numbers of 
SM-like and XM-like quark candidates. 



_ Expected number 
frocess 

of quark candidates 


Correlation 


SM 


XM 


W+HF 


19.5 ±6.4 


0.5 ±0.0 


9.7±3.2 


9.7±3.2 


QCD fakes 


5.4±4.8 


0.48 ±0.06 


2.6±2.3 


2.8±2.5 


Diboson 


2.0±0.4 


0.5 ±0.0 


1.0±0.2 


1.0±0.2 


Mistag 


2.8±0.8 


0.5 ±0.0 


1.4 ±0.4 


1.4±0.4 


Single top 


4.4±0.5 


0.51±0.01 


2.3±0.3 


2.2±0.3 


Total background 


34±8 


0.50±0.01 


17±4 


17±4 


Signal 


700 ±120 


0.562™H 394 ±66 

±0.011(syst) 


306±51 



■ W+HF 

■ Mistag 




Lepton p t (GeV/c) 

FIG. 5: Lepton pj distribution. Shaded histograms show signal and 
background predictions stacked to form the total prediction. 



Monte Carlo modeling of geometrical and kinematic accep- 
tance includes effects of parton distribution functions (PDFs), 
initial- and final-state radiation, and jet energy scale. The PDF 
uncertainty is estimated by varying the independent eigenvec- 
tors of the CTEQ6M H PDF set, varying fheOCD scale 
(Aqcd), and comparing the nominal CTEQ5L 113011 PDF set 
with MRST72 [38]. We vary the parameters that govern 
initial- and final- state radiation and obtain the corresponding 



uncertainty by comparing the results with the nominal one. 
Similarly the uncertainty coming from jet energy scale is esti- 
mated by varying the scale within its uncertainties. An addi- 
tional systematic source comes from the choice of the genera- 
tor (and in particular the hadronization model), for which we 
compare PYTHIA with HERWIG. 

All of these systematic uncertainties affect the predicted 
numbers of signal and background events (for details see Ref. 
j36lo and the efficiency and purity of the pairing and JetQ 
algorithms. An additional systematic uncertainty affects the 
pairing: the choice of the top-quark mass used in the simu- 
lated sample and in the % 2 constraint. We measure this uncer- 
tainty from the shift of the values obtained when comparing 
the nominal results to those from two extra samples generated 
with top-quark masses of 170 and 175 GeV/c 2 . Finally, for 
the JetQ purity systematic uncertainty, we take the value ob- 
tained from the calibration in data and add in quadrature the 
effect of initial- and final-state radiation, since these may be 
different between a bb and a tt environment. In Table ITTT1 we 
show the systematic uncertainties on the pairing efficiency and 
purity, and on the JetQ selection efficiency and purity. These 
systematic uncertainties are assigned only to the signal as for 
backgrounds the statistical uncertainty is dominant. 
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TABLE III: Summary of systematic uncertainties (in %) for the % 2 selection and JetQ efficiencies, and for the pairing and JetQ purities. The 
(0.7) figure is given as information but not used in the total uncertainty since the JetQ purity is calibrated in data and the corresponding scale 
factor already corrects for the Monte Carlo hadronization model. The total uncertainty is calculated as a sum in quadrature of the individual 
uncertainties coming from the different sources. 



Systematic (in %) 



X selection efficiency JetQ efficiency Pairing purity JetQ purity 



Jet energy scale 


0.2 


0.04 


0.1 


0.1 


Initial- and final-state radiation 


0.5 


0.1 


0.2 


0.2 


MC generator 


0.2 


0.1 


0.1 


(0.7) 


Top-quark mass 


0.4 


0.2 


0.9 


0.5 


PDF 


0.7 


0.02 


0.1 


0.02 


Total 


1.0 


0.3 


1.0 


0.6 



IX. SIGNAL PURITY DETERMINATION 

In Table [TT] we show the signal purity that leads to the es- 
timation of the expected numbers of SM-like and XM-like 

I 



quark candidates. The purity is a combination of the pairing 
purity and the JetQ purity as follows: 



Ps — fnb SF„b p„b + ( 1 - fnb SF„ b ) [p wb pjQ SFjq + (1 — pm,) (1 — pjQ SFjq)] , 



(3) 



r 



where f nb is the fraction of signal Monte Carlo events where 
we have misidentified the b jet and SF nb is a scale factor that 
accounts for any difference in the rate of misidentified b jets 
between data and simulation. This is the same scale factor de- 
termined in the measurement of the top-quark-pair production 
cross-section using fo-jet tagging Il36ll . The quantity p nb is the 
probability that a signal event with a misidentified b jet will 
be correctly labeled as SM- or XM-like, pw b is the pairing pu- 
rity for cases where the JetQ was defined, and p jq is the JetQ 
purity for the cases where the pairing criterion was applied. 
These three purities are obtained from simulated events. The 
SFjq is the scale factor between data and Monte Carlo simula- 
tion for the JetQ obtained from the data calibration study (see 
Sec. IVI At . Table iTVl shows the values used in Eq. (0), with 
uncertainties propagated from those in Table IHT1 



is provided in Table [TT] Finally, Table [V] summarizes the im- 
portant analysis inputs to the statistical extraction of results 
described in the next section. 

TABLE V: Estimated numbers of background and signal candidates 
together with the corresponding purities. 



N s 700 ±120 

N b 34 ±8 

p s 0.562 ± 0.004(stat) ±0.011 (syst) 

p b 0.50 ±0.01 



STATISTICAL TREATMENT 



TABLE IV: Inputs to the signal purity. 

7^ 0.079 ±0.001 

SF nb 1.01 ±0.03 

p„b 0.50 ±0.01 

p wh 0.833 ± 0.001 (stat) ± 0.008(syst) 
p JQ 0.608 ±0.001 (stat) ±0.003(syst) 
SFjq 0.99 ±0.01 (stat) ±0.03(syst) 



The equivalent of signal purity for background events is the 
correlation between JetQ and the primary-lepton charge, and 



Once we apply the pairing and JetQ selection to the data, 
we classify each data pair as SM-like or XM-like, and de- 
fine /+ to be the fraction of SM candidates among the data 
pairs. The aim of the measurement is to test the SM hypoth- 
esis (/ + = 1) against the XM hypothesis (/+ = 0). We write 
the likelihood as the product of two Poisson probabilities for 
the observed numbers x + and x~ of SM- and XM-like can- 
didates (respectively), and four Gaussian constraints on the 
nuisance parameters y s , y b , z Ps , and z Pb (the numbers of sig- 
nal and background candidates and the purities of signal and 
background, respectively): 
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FIG. 6: Distribution of the maximum-likelihood estimate of the SM 
fraction / + from pseudoexperiments under the XM (dashed line) and 
the SM (solid line) hypothesis. The arrow shows our result. 



1 w 



10-V 



io- 2 - 



1 3 r ^\ 
10 5 10 4 10 3 10 2 

°-XM 

FIG. 7: Variation of Pxm (the probability of accepting a false XM) 
with clxm ( me probability of rejecting a true XM). The square rep- 
resents our a priori choice of ocxm =1%, corresponding to Pxm = 
0.16%, while the triangle represents the observed p- values and is 
plotted at the coordinates (pxm^Psm)- 



where N+ and N- are the predicted numbers of SM-like and 
XM-like candidates, and N s , Nb, p s , and pb are independent 
estimates of the nuisance parameters (see Table |V). The ex- 
pectations N+ and 7V_ are computed using the following equa- 
tions: 

N+ = Zpjsf+ + {l-z Ps )y s {l-f+)+Zp b yb , (5) 
AL = (1 -z Ps )y s f+ +z Ps y s (l -f+) + {\-z Pb )y b . (6) 



In Fig. [6] we show the distribution of the maximum- 
likelihood estimate / + of / + , as obtained from pseudoexperi- 
ments based on either the SM hypothesis or the XM hypoth- 
esis. We compute two p-values based on / + as test statis- 
tic: psm (pxm) - the probability of observing a value of the 
test statistic as in data or smaller (larger) assuming that the 
SM (XM) hypothesis is true. To reject the SM we require 
Psm < where Usm is the standard 5-sigma discovery 
threshold of 2.87 x 10~ 7 . To exclude the XM we similarly re- 
quire pxm < &XM- We note that increasing (Xxm makes it eas- 
ier to exclude the exotic model, but reduces the exclusion con- 
fidence level 1 — Uxm- To optimize the choice of ocxm while 
taking into account the sensitivity of the measurement, we 
generate pseudoexperiments to compute the probability Pxm 
of not excluding the XM when the SM is true, as a function 
of (Xxm (Fig. [7}. Using this curve we set a^M=l%, slightly 
above the value for which $xm(Mxm)=V-xm ■ 

We also quote a measure of evidence based on the data ac- 
tually observed in the form of a Bayes factor BF, which is the 
ratio of posterior to prior odds in favor of the SM. The BF 
can also be written as the ratio of the likelihood of the SM to 
the likelihood of the XM. The numerator and denominator are 
separately integrated over uniform priors for the nuisance pa- 
rameters. The quantity 21n(BF) can be interpreted according 
to a well-established scale 13911 . 



XI. RESULTS AND DISCUSSION 

In Table [VI] we show the number of events and candidates 
after applying the pairing and JetQ selection and the number 
of candidates corresponding to the SM and XM hypotheses. 



TABLE VI: Observed number of events before and after the pairing 
requirement, observed number of quark candidates with identified jet 
charge, and observed SM-like and XM-like candidates. 



Number of events 


Quark candidates 




Observed After pairing 


JQ denned SM 


XM 


815 397 


774 416 


358 



Candidates whose W-boson charge is opposite to the JetQ 
value are classified as SM candidates, while candidates whose 
W-boson charge has same sign as the JetQ are assigned as 
XM candidates. Figure [8] shows the graphical representation 
of these numbers, where candidates (and SM expectations) 
are distributed as function of the product of the JetQ value 
and the charge of the W boson. Using these numbers we get 
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FIG. 8: Distribution of the product of the W -boson charge times the 
JetQ value. Shaded histograms show signal and background predic- 
tions stacked for the total prediction. SM-like candidates are on the 
negative side of the plot while XM-like candidates are on the posi- 
tive side. The outermost bins correspond to the cases where JetQ is 
exactly ± 1 . 



the profile log-likelihood curve shown in Fig. [9] The mini- 
mum of the curve is at a value of /+ = 0.83. This corresponds 
to a p-value of 13.4% under the SM hypothesis (see red tri- 
angle in Fig. [7j and indicates consistency between CDF data 
and the SM. The p-value under the XM hypothesis is 0.014%, 
which is interpreted as a 99% C.L. exclusion of the XM hy- 
pothesis. The previous measurements have excluded the XM 
hypothesis with at most 95% C.L. [a, 13]. We obtain a value of 
21n(BF) = 19.6, which, according to the interpretive guide- 
lines of Ref. 1331 . constitutes very strong evidence in favor of 
the SM and against the XM. 




FIG. 9: Distribution of twice the negative logarithm of the profile 
likelihood as a function of the fraction of SM candidate events in 
data. 



Table IVHl lists the analysis results for electrons and muons 
separately. In Figs. [lO]and[TT]we show the distribution of f + 
for electrons and muons respectively. Due to its dependence 
on sample size, the measurement sensitivity is lower in each 
lepton subsample than in the full data sample, and an appro- 
priate value of (Xxm is 5% in this case. The XM hypothesis is 



excluded at the 95% C.L. using the electron or muon subsam- 
ple. 

TABLE VII: Number of observed candidates and results for the elec- 
tron and muon candidates separately. 



Electrons 



Muons 



Number of candidates: 

/+ 

PSM 
PXM 

N s 
N b 

Ps 

Pb 



206 SM and 155 XM 210 SM and 203 XM 
1.11 0.57 
65.9% 5.2% 
0.04% 0.7% 
308±51 392±67 
17±5 17±4 

0.56±0.01 0.56±0.01 

0.50±0.02 0.50±0.01 




FIG. 10: Distribution of the maximum-likelihood estimate of the SM 
fraction / + from pseudoexperiments under the XM (dashed line) and 
the SM (solid line) hypothesis for electrons only. The arrow shows 
our result. 

For the muon subsample the p-value under the SM hypoth- 
esis is only 5.2%, compared with 65.9% for the electron sub- 
sample. A % 2 test of the hypothesis that the ratio of XM to SM 
candidates is the same in both subsamples yields a /?-value of 
about 9%, consistent with the discrepancy being a statistical 
fluctuation. 



XII. CONCLUSION 

We present a measurement of the top-quark electric charge 
that relies on the jet-charge algorithm as an estimator of the 
electric charge of high-pr b jets. The measurement uses tt 
pairs reconstructed in final states with one W boson decay- 
ing hadronically and the other leptonically, from a data set 
corresponding to 5.6fb~' of pp collisions collected by the 
CDF II detector. It provides the most sensitive results to sup- 
port or exclude the exotic-quark hypothesis with —4/3 elec- 
tric charge. Our results exclude this hypothesis at 99% C.L.. 
As an additional measure of evidence, the Bayes factor ob- 
tained, 21n(BF) = 19.6, supports very strongly the SM over 
the exotic-quark model hypothesis. 
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FIG. 1 1 : Distribution of the maximum-likelihood estimate of the SM 
fraction / + from pseudoexperiments under the XM (dashed line) and 
the SM (solid line) hypothesis for muons only. The arrow shows our 
result. 
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